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PREFACE 


The  research  described  in  this  report  was  conducted  by  personnel  of  the 
Desert  meteorology  Unit,  of  The  Jacob  Blaustein  Institute  for  Desert  Research, 
Ben-Gurion  University  of  the  Negev,  Sede  Boqer,  Israel  from  1  May  1981  to  30  April 
1982  under  Grant  No.  AFOSR  81-0126  to  the  Ben-Gurion  University  of  the  Negev, 
Research  and  Development  Authority,  P.0.  B.  1025,  Beersheva,  Israel. 

Participating  personnel  concerned  with  the  tasks  described  in  this  report  in¬ 
clude  Prof.  Louis  Berkofsky,  Principal  Investigator,  Dr.  Avraham  Zangvil,  Re¬ 
search  Associate,  Mr.  Zvi  Shipponey,  Mathematician  Programmer,  and  Mr.  Eli  Zangvil, 
Programmer. 

Observational  data  utilized  in  this  study  were  obtained  from  the  Institute's 
micrometeorological  (4m)  tower,  collected  on  a  36-channel  data  logger,  and  ana¬ 
lyzed  in  the  laboratory. 

The  Director  of  the  Institute  during  the  conduct  of  this  study  was  Prof. 

Amos  Richmond. 


This  report  should  be  cited  as  follows: 

Berkofsky, L. ,  1982.  "The  Behavior  of  the  Atmosphere  In  the  Desert  Planetary 
Boundary  Layer,"  Final  Scientific  Report,  prepared  by  the  Desert 
Meteorology  Unit  of  The  Jacob  Blaustein  Institute  for  Desert  Research 
for  the  U.S.  Air  Force  of  Scientific  Research,  AFSC,  Bolling  AFB,  D.C. 
20332. 
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Introduction 

There  exist  a  large  number  of  planetary  boundary  layer  models  each  de¬ 
signed  for  specific  purposes,  (Deardorff,  1974,  Mahrt  and  Lenschow,  1976, 

Stull,  1976,  Heidt,  1977,  Yamada,  1979).  Many  of  these  consider  the  top  of 
the  boundary  layer  a  material  surface.  Some  consider  the  top  of  the  boundary  layer 
to  be  coincident  with  the  inversion,  and  consider  entrainment  across  its  inter¬ 
face.  Some  are  multi-level,  some  are  bulk,  single-level  models.  The  various 
models  are  of  one,  two,  and  three  dimensions.  The  greater  the  number  of  dimen¬ 
sions,  the  greater  the  computational  complexity. 

It  is  possible  to  reduce  the  computational  complexity,  and  still  not 
eliminate  the  three-dimensionality  completely,  by  using  a  variation  of  the 
"momentum  integral"  method  (Schl ichting ,  1968).  By  means  of  this  approach,  the 
vertical  structure  of  several  of  the  variables  is  specified,  and  incorporated 
into  the  vertically-integrated  equations.  In  this  way,  the  model  becomes 
two-dimensional  in  the  horizontal  and  the  vertical  variations  are  incorporated 
in  various  coefficients. 

We  consider  a  model  of  the  planetary  boundary  layer  (depth  approximately 
1  km).  This  layer  is  itself  divided  into  a  surface  layer  (approximately  20  m) 
and  a  transition  layer  which  at  certain  times  and  places  becomes  very  well  mixed, 
and  is  frequently  called  the  "mixed"  layer  (Figure  1).  We  shall  not  insist 
that  this  condition  be  satisfied  everywhere  in  the  model. 

Very  often,  the  top  of  the  planetary  boundary  layer  is  capped  by  an  inver¬ 
sion.  When  this  is  so,  and  when  convection  occurs  below  the  inversion,  the  in¬ 
version  changes  height  due  to  upward  and  corresponding  downward  fluxes  through 
it  by  turbulence.  These  processes  have  to  be  modelled.  Further,  the  processes 
which  we  wish  to  model  are  on  such  a  scale  that  fairly  high  resolution  is  needed-- 
on  the  order  of  10-20  km  in  the  horizontal,  over  a  region  approximately  400-500 
km  on  a  side.  If  then  the  vertical  resolution  is  to  be  very  fine — say  100  m 
(above  the  surface  layer),  the  computation  time  for  solution  of  only  the  boundary 
layer  mesoscale  equations  may  become  prohibitive.  Nevertheless,  the  optimum 
mesoscale  model  must  be  three-dimensional.  As  stated  above,  it  may  be  possible 
to  reduce  the  computation  time  considerably  by  means  of  a  two-dimensional  model, 
and  still  gain  valuable  insights  by  carrying  out  selected  numerical  experiments. 
The  results  of  such  experiments  will  provide  information  for  proper  formulation 
of  future  three-dimensional  mesoscale  models.  For  this  reason,  we  propose  to 
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parameterize  quantities  in  the  vertical  by  assuming  that  we  are  permitted  to 
specify  their  variation  with  height.  If  we  then  integrate  the  equations  verti¬ 
cally  up  to  the  inversion  height,  and  introduce  the  parameterized  variables,  the 
vertical  variation  goes  into  the  coefficients,  while  the  equations  themselves 
apply  to  parameters  at  one  or  two  specific  levels.  The  procedure  will  be  des¬ 
cribed  in  greater  detail  below. 

Model  Equations 

We  shall  concern  ourselves  with  a  form  of  the  primitive  equations  derived 
by  averaging  over  a  horizontal  area  A  x  Ay  which  is  large  enough  to  contain  the 
sub-grid  scale  phenomena,  but  small  enough  to  be  a  fraction  of  the  mesoscale 
system.  We  define 


1 


")  dy.  d 


Cy  r:  Cf  -f~  cy  ' 


(1) 


(2) 


where  ‘Y'is  any  scalar  variable. 

With  the  above  definitions,  the  appropriate  equations  are,  approximately 
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In  the  above  set  of  equations,  the  unbarred  variables  u,  v,  w,  u  ,  v  ,  H,  B, 

9  9 

q  are  all  mean  values  according  to  Equation  (1).  We  have  used  the  Boussinesq 
Assumption,  and  the  variables  are  defined  as  follows: 

u,  v,  w  =  components  of  the  wind  vector 

u  ,  v  =  geostrophic  wind  components 
y  y 

H  =  CpT+Lq+gz  =  moist  static  energy 
q  =  specific  humidity 
Q|p  =  divergence  of  net  radiation  flux  = 
x,  y,  z,  t  =  space  variables  and  time 
Equations  (3)  and  (4)  are  the  horizontal  equations  of  motion.  Equation  (5)  is 
the  continuity  equation.  Equation  (6)  is  the  thermodynamic  energy  equation. 
Equation  (7)  is  the  specific  humidity  equation.  The  lower  boundary  condition  is 

\xr  =  =  terrain  height 

T  •  (8) 

=  given  K  1 

In  the  above  system,  we  have  chosen  to  treat  the  moist  static  energy.  The 
solution  for  this  quantity,  together  with  that  for  the  specific  humidity,  will 
enable  us  to  recover  the  temperature  when  needed.  In  the  present  investigation, 
we  assume  that  any  condensed  moisture  stays  in  the  air.  Thus,  we  do  not  treat 
clouds  or  precipitation  explicitly  in  this  model  (but  implicit  predictions  are 
possible).  We  indicated  earlier  that  changes  of  inversion  height  will  have  to 
be  modelled.  This  is  particularly  important  in  desertification  studies,  be¬ 
cause  many  desert  areas  are  located  near  the  descending  branch  of  the  Hadley 
cell,  where  inversions  exist.  This  can  be  seen  from  Israeli  data  'Shaia  and 
Jaffe,  1976),  where  there  are  222  days,  on  the  average,  per  year  of  mid-day 

-4- 


j 


inversions  just  about  100  km  north  of  the  beginning  of  the  Negev  desert. 

In  order  to  expedite  deriving  appropriate  expressions  for  predicting 
inversion  height,  we  first  derive  inversion  "interface"  conditions.  The  re¬ 
sults  are  essentially  the  upper  boundary  conditions.  Let  h(x,y,t)  be  the  in¬ 
version  height.  Let  ^  be  a  small  layer  of  constant  thickness  above  the  in¬ 
version  level.  In  Mahrt  and  Lenschow  (1976)  this  is  called  the  turbulent  in¬ 
version  or  transition  layer,  which  is  sufficiently  thin  so  that  terms  of  0  (  l  ) 
may  be  neglected  in  the  integrated  equations. 

(  r-*  +  £ 

Define  (  7  =  —  j  (  7  d  X  (9) 

and  assume  that  w  may  be  replaced  by  its  value  uu^at  the  bottom  of  layer. 


Let  A  (  ~)  =  (  -  (  \ 

By  Leibniz's  rule. 


^  , _ 


where  "?(-•  =(>,7,0 

We  apply  Equation  (9)  to  Equations(3)  -  (6)  inclusive,  and  obtain 


(10) 


(11) 


(  lA 

V. 


A  U  4-  A 


(12) 


LO-^  )AV  f  4  (va-O-7 
/  ' 


=  -  (Vu>:\ 


(13) 


Each  of  the  above  equations  can  be  viewed  as  parameter izations  of  the  ver¬ 
tical  eddy  fluxes,  or  as  prediction  equations  for  h  if  these  eddy  fluxes  are 
Known.  The  quantity  -  wh),  which  represents  the  motion  of  the  air  relative 

to  the  inversion,  is  called  the  "entrainment  velocity,"  w  .  w,  is  the  larqer- 

e  h  .  vo 

scale  vertical  velocity  at  the  interface.  The  terms  involvinq  and  -~x 

aT-  d  y 

are  usually  omitted  in  derivations  of  these  interface  conditions,  since  most 


inversion  height  models  assume  horizontal  homogeneity. 

Our  approach  will  be  to  integrate  the  system  of  Equations  (3)  -  (7)  in¬ 
clusive,  with  a  modified  treatment  of  Equation  (5),  with  respect  to  z  from  z  =  k 
(constant)  =  height  of  surface  layer,  to  z  =  h  (inversion  height),  and  then  to 
introduce  modelling  assumptions  for  the  variables,  i.e.,  to  specify  their  vari¬ 
ations  with  height.  If  we  do  this,  it  becomes  possible  to  express  all  of  the 
jump  quantities  in  Equations  (12)  -  (15)  inclusive  in  terms  of  their  values  at 
a  specific  level . 

Parameterizations 

For  the  winds  in  the  surface  layer,  we  assume  a  constant  flux  profile 


u(?)  r 


(16) 


where  k1  =  von  Karman  constant  =  0.4,  U*=  friction  velocity,  zq  =  roughness 

parameter . 

In  the  transition  layer,  we  assume 


t)  =  BW  v  (i, y,  t') 


(17) 
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where 


A 


(  j  a  t 


(18) 


We  derive  a  modelling  approximation  for  temperature  over  land  above  which 
a  nighttime  inversion  exists,  and  where  surface  heating  occurs  during  the  day, 
leading  to  convection  and  turbulent  mixing.  In  Figure  1,  an  inversion  exists 
in  the  early  morning.  The  temperature  at  z  =  k  is  TeLT  where  L  means  land  and 
I  means  initial  value.  The  temperature  increases  linearly  with  height  according 
to 


TGy,  0s)  -  T^li  -t-  T,  (19) 

is  a  constant  lapse  rate. 

As  convection  sets  in  and  h  increases,  the  temperature  changes  with  height 
according  to 


ht'j)  "2 1  t)  -  Tfcu  ('t't'ji  0  lA  (20) 

where,  the  lapse  rate  Y"(z)  need  not  be  constant. 

Thus,  the  temperature  jump  is  seen  to  be 

^  T  ~  (T^x  -  -f  ^  ^  ( A ~§k  )  4-  T  £  (21) 

Suppose  that  air  flows  from  land  to  water,  that  the  water  is  cooler  than  the 

land,  and  that  the  air  above  the  water  is  stably  stratified.  As  the  warm  air 
moves  over  the  water,  it  is  cooled  from  below,  i.e.,  there  is  a  downward 
transfer  of  heat  from  the  warm  air  to  the  cooler  water.  At  the  same  time,  tur¬ 
bulent  mixing  is  taking  place,  tending  to  destabilize  the  temperature  profile. 
After  the  mixing  process  sets  in  the  temperature  changes  with  height  according 
to 


’■  ^  ?  ^  -  T4  ^ ,  t)  +  (2  )  [-£  K  ^  t)  -  4 

The  inversion  strength  (temperature  jump)  over  water  is 

4T  ^  (l^L1  -T*  j  +  -  -rtt)  ]  7-4)  +  V;  S> 

Similarly,  the  mixing  ratio  profiles  are  given  by 

C'r  (>  i*y,  4-  S' 1  (  2  -"H) 

where  &  is  constant  =  initial  lapse  rate  of  moisture 

^■7'^  -  7,7,7  +  5(v>  7-4) 

where  ^(2)  =  lapse  rate  of  moisture  over  land 

*r  (*,'/, l, t)  ~  fay,  b)  f 

wliere  lapse  rate  of  moisture  over  water. 

The  moisture  jump  is  then 


Then  the  moist  static  energy  is 
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and  the  moist  static  energy  jump  is 


(t  -l)  Wc^+lOS 


(29) 


We  shall  apply  the  operator  Equation  (18)  to  the  basic  Equations  (3)  -  (7) 
inclusive,  introduce  the  modelling  assumptions  Equations  (17),  (19)  -  (24)  in¬ 
clusive,  and  use  the  interface  conditions  Equations  (12)  -  (15).  However , 
these  will  only  be  used  to  eliminate  (u'uy1^  and  ("v1  ur1),^  from  the  first 
two  equations  of  motion  Equations  (3)  and  (4).  The  first  law  of  thermodynamics 
becomes  very  complex  if  we  use  that  procedure  there.  Instead,  we  introduce  the 
closure  assumption 


(  tAT  H,N)£  -  (30) 

i.e.,  the  turbulent  flux  of  heat  at  interface  height  is  a  fraction  of  that  at 
the  top  of  the  surface  layer,  and  oppositely  directed.  This  assumption  has 
been  widely  used  in  the  literature  (see  Carson  (1973)).  We  shall,  in  fact  use 
the  same  closure  in  the  interface  condition  Equation  (14)  for  AH.  The  resulting 
equation  will  then  become  the  prediction  equation  for  h. 

For  the  fluxes  at  2  =  we  write 

W  UJ')^  ~  ~  ^ 

|v*\ 
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PW  =  P^cyC^  [ft  (A)  u](T^-T5)  (33) 


In  the  above,  C*  is  a  turbulent  transfer  coefficient  for  momentum,  CHc,  is 
a  turbulent  transfer  coefficient  for  heat,  T^is  ground  temperature,  V/^  = 
ground  soil  moisture.W^  potential  saturation  value  of  V/ 
lie  shall  also  parameterize  by 

fc 


0 


f 


S  4  rT-C  [ft  fa,  *  b.e  Q  _  ,  +lss. 

_ _ _ - _ Tft  , 

~ . e*  ft-t) 


«  (35) 

_  [mveT 

(j/'K-'fi) 


Here  ^  -V i s  incoming  solar  radiation,  C^is  ground  emissivity,  a^  and  b^  are  con¬ 
stants  associated  with  downcoming  infrared  radiation,  e  is  saturation  vapor 
pressure,  is  ground  albedo,  defined  by 


R**  --  *  +  b  ,  b  Co 

a,  b  constants  | 


(36) 
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Thus  ^decreases  as  soil  moisture  increases. 

It  is  seen  from  the  above  that  we  shall  need  equations  for^c^and  V- 
These  will  be  adapted  from  Deardorff  (1978). 

By  applying  all  of  the  modelling  assumptions,  we  obtain  the  following: 
Prognostic  Equations 

These  are  written  for  land,  for  water,  we  replace  X  by  of ,  £  by  £  • 
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First  Equation  of  Motion 
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-  ~  Q 


U-4) 


Second  Equation  of  Motion 
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First  Law  of  Thermodynamics 
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Moisture  Equation 
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Inversion  Height  Equation 


where 

a2  -  ^[r(-rc^] 

a2  =  Pt(A+s) (cyr, iU,)- M) \Sf  *(*)  f 
'  a3  *  BU+S)  (CyV.tLO  -S^^ra)  +iR^] 

a4  -  \.r  v\-fe.n "  ^0  +  c^rj  +lO0 

a5  ~  (-^+0  U-feo:  ~fV{'£)  -\-  A'Cr^'h^  y  ^  -f-  L  O  £ 

aG  =  i&  U  -ft')  H-ftx  “  &  PC  4-  6  (^+f)(cp  H  +  L-S,')  J 
In  the  above  derivation,  we  have  neglected  aF. 
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Ground  Temperature  Equation 


(42) 


where 


i 


■ir  __ 

ha  =  +  L 

p  -  soil  density  -r  „  . , 

Is  'i.  -  deep  soil  temperature 

C5  =  soil  specific  heat 
<i  i  * 

$5  .=  soil  thermal  diffusivity 
^  -  period  of  1  day 

Equation  (37)  is  referred  to  as  the  "force  restore"  equation  for  ground  temp¬ 
erature 


cv  ( -  vO 
P-  t 


(43) 


Soil  Moisture  Equation 

where  and  C2  are  constants 

(?w  =  density  of  liquid  water 
P  =  precipitation  rate  (prescribed) 
bulk  moisture  (analogous  to  T*.) . 

Notice  that  all  capped  coefficients,  i.e.,  A  etc.,  are  functions  of  h. 

We  still  need  diagnostic  equations. 

Diagnostic  Equations 

The  parameterizations  themselves  are  diagnostic  equations.  We  have  already 
stated  an  empirical  diagnostic  equation  for  ground  albedo  (Equation  (31)).  We 
still  need  expressions  for  the  vertical  velocity  w„  and  the  geostrophic  wind 
values  Ug-,  xr^  . 

When  deriving  the  appropriate  continuity  equation  from  Equation  (5)  and  the 
lower  boundary  condition  Equation  (8),  we  have  to  distinguish  between  two 
cases  -  when  tT  <==K  ,  and  1=r  >\  .  Thus,  vertical  integration  of  the  continuity 
equation  and  the  application  of  Equations  (16)  and  (17)  leads  to 
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Equations  of  Continuity 


can  be  obtained  by  setting  It  =  h  in  the  upper  limit  of  the  integrals  from 
the  appropriate  equation. 

In  terms  of  the  modelling  assumptions,  we  have 


f  ^  Ulx  +  (i) 

±(Du-o~ 


Mean  Geostrophic  Wind  Components 


The  system  of  prognostic  Equations  (37)  -  (43)  inclusive,  together  with 
the  diagnostic  Equations  (33),  (34),  (35),  (44),  (45),  (46),  comprise  the  forecast 
scheme  for  prognostic  variables  u  ,  \r/  and  the 

diagnostic  variables  (  (  (  ,  E^_  ,  provided  t  i  and 

the  various  functions  of  1  are  known. 

It  is  possible  to  add  a  dust  concentration  equation  to  this  system,  for 
example,  a  time-dependent  equation  in  which  the  time  rate  of  change  of  dust 
concentration  is  given  by  a  balance  (or  imbalance)  among  sedimentation, 
diffusion,  and  transport. 

A  hierarchy  of  models  based  on  the  above  general  system  can  be  envisioned. 

A  Hierarchy  of  Models 

The  simplest  time-dependent  model  which  may  be  envisioned  from  the  above 
general  model  is  one  in  which  there  is  no  horizontal  advection,  and  hence  no 
horizontal  divergence  and  no  vertical  motion.  We  consider  a  dry  model  over 
flat  terrain. 

^  +-  ^  =-C4^A(-fO]  lv|u  (47) 

^  (k-'k) 

1st  Equation  of  Motion 


\ir 

'J'f 


+  -f  (u- 4  [i -  icbjll  ~  jb --  -g,  I BtfTl  V 1 V  / 

p  To  T\  * 


(48) 


2nd  Equation  of  Motion 


>t>  fc,-r s)  (49) 

>tr  P  ('R-'R') 


First  Law  of  Thermodynamics 
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X(4-4.)  =  1  A^u  (Ttlt-Ty)  ) 


(50) 


Inversion  Height  Equation 


HI Cft 


vA 


P>M, 


L 


(51) 


Ground  Temperature  Equation 

Notice  that,  if  h  =  const.,  the  two  equations  of  motion  reduce  to  the 
E<man  equations.  According  to  Equation  (46)  neglect  of  horizontal  gradient 
leads  tot^s'its  o.  However,  we  should  include  a  basic  current,  and  this  is 
done  by  writing 


U|  -  constant 


(52) 


The  next  set  in  the  hierarchy  will  be  derived  by  adding  horizontal  trans¬ 
port  in  the  x-direction  and  vertical  motion  to  the  above. 
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1st  Equation  of  Motion 


(53) 
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2nd  Equation  of  Motion 
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First  Law  of  Thermodynamics 
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Inversion  Height  Equation 


~  0 


(56) 


(57) 


Equations  of  Continuity 
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Ground  Temperature  Equation 

A 

According  to  Equation  (46),  =  0  unless  we  have  a  constant  basic  flow. 


i  .e. 


/V 


>Uc 


U£  ~ 


A.  1T5 

-f  *1 


_i2_  '^ur  =  Cov\st$Y\t 


A,  -A  A-x  C-i-A) 

L- 


(59) 


A.  A 

Equations  (53)  -  (59)  inclusive  solve  the  system  for  u  ,  'V  ,  T^(  -k  , 

Tea.  ,  WX  ,  .  £  ,  A-  ,  A-  K+ft  ^  ,  r  (-O 

are  all  functions  of  h-, 

■fi ,  ?,  ^  ,  Cjl  ,  e\ 

are  all  constants,  STis  given.  The  functions  A(5),  B(‘?),f(^)  are  empiri¬ 
cally  determined. 

These  are  the  next  systems  in  the  hierarchy  of  those  we  plan  to  investi¬ 
gate.  Following  these  we  will  incorporate  the  moisture  in  one-dimension,  and 
the  terrain  in  one  dimension.  We  will  then  go  on  to  two-dimensional  (space) 
models . 

A  Test  of  the  System 

A  very  simple  test  of  the  system  can  be  made  if  we  consider  an  even 
simpler  version  of  the  inversion  height  equation.  If,  in  Equation  (56)  we 
omit  horizontal  advection  effects,  we  obtain 
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I  C»«,  A  ft')  0  (Tgc  -T«)| 
'  [iT-H#)]  [■&-■£) 4(Tir^T^ 


(60) 


This  equation  is  similar  to  that  of  Tennekes  (1973). 

Of  course,  to  test  Equation  (60),  prediction  equations  are  required  for 

/s 

^GR’  Tk’  ^ »  anc*  a  diagnostic  equation  for  What  we  have  done  is  to  use 
real  data  for  the  heat  flux  Lt(TGR  -  T^),  and  the  temperature  anomaly 

(Tki  -Tk).  For  ,  we  have  assumed  values  shown  by  the  curve  in  Figure  2. 

The  data  are  for  Sede  Boqer,  Israel  for  18  June  1979.  We  have  assumed  an 
initial  value  of  h  =  300  m  at  0900  LST.  We  use  =  0.2,  =  0.5  x  10”4  deg 

cm-1,  y(h)  =  -0.5  x  10-4  deg  cm'1,  C^q  =  1.695  x  10"3,  k  =  20  m. 

The  results  are  shown  in  Figure  3.  On  the  same  figure  is  the  march  of  heat 
flux  near  the  ground.  Although  no  verification  data  are  as  yet  available,  it 
can  be  seen  that  the  variation  of  inversion  height  throughout  the  day  is  rea¬ 
sonable.  The  maximum  of  1051  m  is  reached  at  1800.  This  maximum  lags  the  heat 
flux  by  about  6  hours,  which  is  in  agreement  with  the  O'Neill,  Nebraska  re¬ 
sults  shown  in  Figure  1  of  Carson  (1973).  Our  results  are  very  sensitive  to 
vertical  velocity  variations.  Figure  3  shows  the  evolution  of  the  inversion 
height,  using  the  same  data  as  above,  but  with  wh  =  0,  (see  Equation  (60)). 

It  is  clear  that  the  larger-scale  vertical  velocity  plays  an  important  role  in 
determining  inversion  height,  especially  at  night. 

In  the  near  future,  we  plan  to  release  mini-sondes  in  order  to  obtain 
verification  data.  We  will  also  attempt  to  obtain  data  by  means  of  a  tethered 
balloon.  We  have  an  AN/TPQ-11  vertically  pointing  radar,  which  is  capable  of 
detecting  inversions.  We  hope  to  have  this  instrument  operable  during  the 
coming  year. 

We  have  also  carried  out  calculations,  based  on  Equation  (60),  to  test  the 
effect  of  a  water  body  of  the  air  flowing  over  it.  In  this  case,  we  must  of 
course  include  moisture.  Using  the  various  modelling  approximations,  the  in¬ 
version  height  equation  becomes  (recall  that  WgR  =  W^  over  water), 
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(L r'  -  UJC>\  +  +jfr*y +-Wpuy) 


Here,  if  we  wish  to  feed  in  data  over  water,  we  are  in  a  worse  position  than 
previously,  as  we  have  no  such  data  for  the  Sede  Boqer  area.  We  must  there¬ 
fore  make  some  heavy  assumptions.  Suppose  we  assume  that 


T^-T^)  +-  =fc*-T*),  («> 


[(y-y  +  =  (t,x- 


This  is  based  on  the  fact  that  both  (T^  -  T^)  and(T^  -  T^)  are  smaller  than 
their  values  over  land.  However,  the  latent  heat  term  increases  their  values. 
Then  the  main  difference  between  land  and  water  lies  in  the  factor  Cr,  —  9-CfOj 
♦  y  t? ,  -  &  cy} .  D.  We  can  generate  a  series  of  results  for  various  val- 


% 
ues  of 


a,  .ia,  .... ,  z.O'  ^ \  0-1- 
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Of  course,  if  we  keepW^  the  same  as  that  over  land,  we  are  assuming  the  same 
kind  of  convection  over  water  as  over  land.  The  calculations  can  also  be 
carried  out  for  various  w^. 

We  have  carried  out  one  calculation,  with  w.  as  given  in  Figure  2.  We 
-4  n 

have  used  D  =  1.2  x  10  .  All  other  conditions  were  the  same  as  previously 
The  resulting  curve  for  h  is  shown  by  the  middle  -|-jne  -jn  Figure  3.  The  entire 
curve,  after  the  first  two  hours,  lies  below  that  for  land.  The  maximum  of 
950  m  is  reached  at  1800,  as  before. 

The  principal  conclusion  to  be  drawn  is  that  the  inversion  undergoes  an 
evolution  over  water  similar  to  that  over  land,  but  that  it  reaches  a  lower 
maximum  than  that  over  land.  Furthermore,  the  minimum,  in  the  early  morning 
hours,  was  also  lower  than  over  land.  This  result  has  important  consequences 
for  problems  like  radar  ducting,  dust  spreading,  and  air  pollution. 

In  order  to  test  this  result  further,  it  will  be  necessary  to  obtain  data 
from  a  meso-net,  both  over  land  and  over  water.  Then  the  system  of  Equations 
(53)  -  (59)  inclusive  can  be  integrated  for  a  thorough  test. 

Conclusions 

We  have  derived  a  general  system  of  vertically  integrated  equations  for  the 
desert  planetary  boundary  layer.  From  this  system,  a  hierarchy  of  models  can 
be  designated.  We  have  tested  the  simplest  possible  version  in  the  hierarchy, 
where  real  desert  data  were  fed  in  to  the  inversion  height  equation.  The 
evolutionof  the  desert  inversion  was  very  realistic,  although  the  results  have 
not  yet  been  verified.  We  have  shown  that  the  vertical  velocity  of  the  meso- 
scale  motions  at  inversion  height  is  an  important  parameter  in  determining  in¬ 
version  oscillations.  This  is  particularly  true  at  night,  when  sinking  motion 
tends  to  depress  the  inversion.  The  effect  of  vertical  motion  is  clearly 
brought  out  when  the  calculations  are  carried  out  for  zero  vertical  motion. 

Then,  the  large  upward  flux  of  heat  raises  the  inversion  during  the  day,  but 
the  downward  flux  during  the  night  is  insufficient  to  lower  the  inversion  height, 
which  remains  constant  after  reaching  a  maximum  in  the  late  afternoon.  Using 
assumed  data  for  water,  the  .inversion  height  simulation  appears  realistic,  and 
more  or  less  parallels  that  over  land.  But  the  maximum  height  is  less  than  that 
over  land,  as  is  the  minimum. 
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Recommendations 

It  is  reconmended  that  we  continue  with  the  next  set  in  the  hierarchy. 

At  the  same  time,  we  plan  to  release  mini-sondes,  launch  a  tethered  balloon, 
and  put  into  operation  a  vertically  pointing  radar,  which  is  capable  of  de¬ 
tecting  inversions.  We  also  plan  to  gather  data  from  a  meso-net  for  additional 
verification. 

We  are  purchasing  an  instrument  for  sampling  the  dust  content  of  the 
desert  atmosphere.  This  instrument  will  separate  the  dust  into  five  sizes. 

We  will  be  able  to  obtain  the  individual  and  total  size  concentration.  These 
data  will  serve  as  input  into  a  numerical  dust  model,  and  will  also  serve  as 
ground  truth  information  for  lidar  measurements  giving  dust  distribution  as  a 
function  of  height.  Details  of  the  derivation  of  a  dust  model  equation  and 
some  preliminary  results  are  given  in  a  Proposal  to  the  Department  of  the 
Air  Force  (AFSC)  to  extend  subject  grant. 
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Figure  2.  Assumed  values  of 


Figure  3.  Inversion  Heights  (h)  and  Surface  Heat  Flux  (HF)  vs.  Time 


Heat  Flux  (cm-s_l-deg  x  10  2) 


3 

CD 


zr  — 


ro 

OJ 

-f* 

cn 

cr> 

~n! 

00 

<D 

O 

o 

o 

O 

o 

o 

o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

Inversion  Height,  h(m) 


h(land) 


